Dependence of the transition from Townsend to glow discharge on secondary emission 
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In a recent paper Sijacic and Ebert have systematically studied the transition from Townsend to 
glow discharge, refering to older work from von Engel (1934) up to Raizer (1991), and they stated 
a strong dependendence on secondary emission 7 from the cathode. We here show that the earlier 
results of von Engel and Raizer on the small current expansion about the Townsend limit actually 
are the limit of small 7 of the new expression; and that for larger 7 the old and the new results vary 
by no more than a factor of 2. We discuss the 7-dependence of the transition which is rather strong 
for short gaps. 



In the recent article [i| , the transition from Townsend 
to glow discharge was re-investigated with analytical and 
numerical means. On the analytical side, a systematic 
small current expansion about the Townsend limit was 
performed and it was stated: 

"The result agrees qualitatively with the one given by 
Raizer and Engel and Steenbeck In particular, 
the leading order correction is also of order a"{j/^Y. 
However, the explicit coefficient of p differs: while the 
coefficient in 2^ _3] does not depend on 7 at all, we find 
that the dependence on 7 is essential, as the plot of F 
in Fig. 1 (of 01) clearly indicates. In fact, within the 
relevant range of 10"^ < 7 < 10°, this coefficient varies 
by almost four orders of magnitude. We remark that it 
indeed would be quite a surprising mathematical result 
if the Townsend limit itself would depend on 7, but the 
small current expansion about it would not." 

Here we remark that while the new systematic calcu- 
lation was correct, the interpretation and comparison to 
earlier work requires some correction. 

To be precise, the model treated in 0, 0, and by 
many other authors is a one-dimensional time indepen- 
dent Townsend or glow discharge characterized by the 
classical equations 

d^Je - I Jel a(l^l), d^J+ = I Jel a{\E\) , (1) 



dxE = — (n+ - ne), 
for electron and ion particle current 



(2) 



(3) 



and electric field E. Impact ionization in the bulk of the 
discharge is given by the Townsend approximation 
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(In 0, the generahzed case Q;(|-E|) = e-xp {— Eq / \E\Y 
was treated.) Since ions are generated by impact within 
the gap and drift towards the cathode, there are no ions 
at the anode at a; = 



The ions impacting on the cathode liberate free electrons 
with rate 7, therefore secondary emission from the cath- 
ode is given by 



\Je{d)\^-l\J+{d)\, 



(6) 



where a; = d is the position at the cathode. The electrical 
potential difference between the electrodes is 
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dx E{x), 



and the total electric current is 



J = e (n^ 



Ueflc) E, dxJ = 0. 
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It is useful to introduce dimensionless voltage and cur- 
rent 
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where j = j / fi with the definition of j from Jj. It should 
be noted that only bulk gas parameters have been used 
as units; therefore the dimensionless u and j are inde- 
pendent of 7. 

Further dimensional analysis yields that the current- 
voltage-characteristics u — u(j) can depend on three pa- 
rameters only, namely on the dimensionless gap length 
L = aod, on the coefficient 7 of secondary emission and 
on the mobility ratio /i = /x+//Xe- In practice, the depen- 
dence on the small parameter n is almost negligibly weak 
therefore u = u{j,L,j). Here the dimensionless gap 
length L is related to pd through L ~ Apd a.s long as the 
coefficient ag is related to pressure like ao = Ap. 

How strongly does the characteristics u = u{j,L,j) 
depend on 7? In 0, Sijacic and Ebert calculated the 
whole Townsend-to-glow regime numerically and derived 
by expanding systematically in powers of current j about 
the Townsend limit 
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which gave an excehent fit to the numerical solutions. 
Here 
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and £t and are field and potential in the Townsend 
limit of "vanishing" current, i.e., with breakdown values 



1 



L 



HL/L,y 



ut 



ln{L/L^ 



(14) 



The minimal potential ut is L^e^, it is attained for ga; 
length L = L^e^ on the Paschen curve ut — ut{L) \ 

In m, it was argued that the coefficient Ase in H1U|) 
strongly depends on 7 due to the factor F{'y, n) in (|ll|l . 
This factor strongly depends on 7, for small 7 actually 
in leading order like L^/12. (Note that there is a dis- 
crepancy between equation (50) in for F{'-f,fj,) which 
is reproduced as Eq. ifT^ in the present paper, and the 
plot in Fig. 1 of Q for 10^^ < 7 < 10°. Equation (50) in 
[l| is correct and the figure erroneous. -F(7, /i) actually 
varies by five orders of magnitude on 10^^ < 7 < 10°, 
not only by four.) 

At this point, the question how the remaining factor 
depends on 7 was omitted. In fact, the denominator 
(aSr)^ in H10|l has in leading order the same strong de- 
pendence on 7, since 
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according to the Townsend breakdown criterion aL — 
L-y, cf. (fT^ - flUl . Therefore the leading order depen- 
dence on L"^ of the coefficient of in HlOfl is cancelled 
and replaced by a dependence on , while the term with 
a" has the classical explicit form 

£t a" _1- 2£t _ \n{L/L-,) - 2 
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In Q , another small current expansion was derived 
from assuming n+ ^ Ue and n+(x) w const. 

This approximation was criticised in P|, since it is in 
contradiction with the boundary condition ^ — how- 
ever, for very small 7, it is a good approximation in a 
large part of the gap. The resulting equations (8.8) and 
(8.10) from 2] read in the notation of the present paper 



U 



1 - 2£t 



48 2£t 



J_ 

12 



(17) 



Jl 



2d3 



(18) 



(Here a misprint in was corrected, namely the missing 
factor Ut in the coefficient of in (|17l) . and the factor 
l/(87r) in (8.8) is substituted by eo/2 in (fH^ . since we 
here write the Poisson equation Q in MKS units rather 
than in Gaussian units, cf. (8.6) in 

In l(T7|) . the physical current density J is compared to 
Jl- Jl is the current density where deviations from the 
Townsend limit through space charges start to occur; it 
exphcitly depends on 7 through Ut ^M- 

Comparison of the results of Sijacic/Ebert (SE) 110|) 
and Engel/Raizer (ER) (|17|l show that the coefficients 
AsE.EB^ in the expansion IjlOl) are related like 
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The coefficients Ase and Aer depend in the same way on 
L, and they are essentially independent of /i for realistic 
values of fi. Therefore the ratio Ase/ Aer depends only 
on 7 as shown in Fig. 1. For 7 — > 0, the ratio tends to 
unity. For a large range of 7 values, the deviation is not 
too large, approaching a factor 0.44 for 7 = 10^^. 




FIG. 1: The ratio Ase/Aer of tlie small current expansions 
by Sijacic/Ebert and Engel/Raizer as a function of 7. 



Fig. 2 shows that the factor Ase indeed strongly de- 
pends on 7 for the given L. 

The strong dependence of Ase or Aer. on 7 for a 
given short gap length L means that we can obtain both 
negative and positive differential resistance dU/dJ close 
to the Townsend limit for the same gap length. There- 
fore the choice of 7 is important since it can change the 
differential conductivity and therefore the stability of a 
Townsend discharge in a short gap. 
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FIG. 2: The normalized coefficient A = 24 Ase/{L^ In" L) as 
a function of 7 for gap lengths L — Apd = 15, 30, 60, 120, 
240 (dashed and solid lines with labels). 
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